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Abstract 

In this paper, we propose a general approach for describing hadronic B decays. Using 
this method, all amplitudes for such decays can be expressed in terms of contrac- 
tions, though the matrix elements are not evaluated. Many years ago, Buras and 
Silvestrini proposed a similar approach. However, our technique goes beyond theirs 
in several ways. First, we include recent theoretical and experimental developments 
which indicate which contractions are negligible, and which are expected to be smaller 
than others. Second, we show that all S-decay diagrams can be simply expressed in 
terms of contractions. This constitutes a formal proof that the diagrammatic method 
is rigourous. Third, we show that one reproduces the relations between tree and 
electroweak-penguin diagrams described by Neubert and Rosner, and by Gronau, Pirjol 
and Yan. Fourth, although the previous results hold to all orders in a^, we show that 
it is also possible to work order-by-order in this approach. In this way it is possible to 
make a connection with the matrix-element evaluation methods of QCD factorization 
(QCDfac) and perturbative QCD (pQCD). Finally, using the contractions approach, 
we re-evaluate the question of whether there is a irK puzzle." At 0(a°), we find 

that the diagram ratio \C' /T'\ is about 0.17, a factor of 10 too small to explain all the 
B — * nK data. Both QCDfac and pQCD find that, at 0{al), the value of \C'/T'\ may 
be raised to only about 2-3 times its lowest-order value. We therefore conclude that, 
assuming the effect is not a statistical fluctuation, it is likely that the value of \C' /T'\ 
is similar to its 0{a^) result, and that there really is a i? ^ nK puzzle. 
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1 Introduction 

The study of hadronic B decays'^ has a long and storied history. Much work has gone 
into all aspects of this problem. 

One important development was the operator-product expansion (OPE). Here the 
effective Hamiltonian for quark-level b decays takes the form |lj 

C ( \ 

^e,, = ^ E E A('')(Ci(/i)o?(/i) + C2(/i)of(/i))-A;'^)E^«(/^)o.(/i) , (1) 

V q=d,s \p=u,c i=3 / 

where A^''^ = VpbV*g. /i is the renormalization point, typically taken to be 0{mh). The 
Wilson coefficients Ci (WC's) include gluons (QCD corrections) whose energy is above 
/i (short distance), while the operators Oi include QCD corrections of energy less than 
fi (long distance). All physical quantities must be independent of /i. Note: factors of 
Gpl\pl are omitted for the remainder of this paper. 
The operators take the following form: 

0\ = iPaba)v-A {qpPp)v-A , Of = {pabf3)v~A {qpPa)v-A , (2) 

summed over colour indices a and (3. These are the usual (tree-level) current-current 
operators induced by ly-boson exchange. 

O-i = {qaba)v~A E iqWl3)v-A , O4 = {qJ)p)v-A E {q^V-A , 
q' q' 

O5 = {qaba)v-A E iq'fsq'p)v+A , Oe = {qJ)p)v-A E %q'a)v+A , (3) 

q' q' 

summed over the light flavors q' = u, d, s, c. These are referred to as QCD penguin 
operators. 

3 3 

O7 = 2 iqMv~A E ^9' (qMv+A , Os = - {qabp)v-A E v {ma)v+A , 

q' q' 

3 3 

O9 = 2 {qMv~A E ^1' iq'l5q'l3)v-A , Oio = - {qabp)v-A E ^1' ima)v-A , (4) 

q' q' 

with Cg/ denoting the electric charges of the quarks. These are called electroweak- 
penguin (EWP) operators. The quark current {qiq2)v±A denotes gi7^(l ± 75)g'2- AH 
quark-level b decays can be expressed in terms of the effective Hamiltonian. 

An alternative description is the diagrammatic approach j2], in which quark- level 
b decays are characterized by various topologies. These include the colour-favored and 
colour-suppressed tree amphtudes T and C, the gluonic penguin amplitude P, the 

*In this paper, we present results for B mesons, which contain a b quark. We will nevertheless 
continue to refer to B decays. 
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colour-favored and colour-suppressed EWP amplitudes Pew and -P^w the exchange 
diagram E, the annihilation amplitude A, and the penguin-annihilation diagram PA. 
This formalism has been extensively used, but the relation to the OPE and the effective 
Hamiltonian has not been made entirely clear. 

Of course, the true decays involve mesons, and do not take place at the quark 
level. That is, the decay is i? — > F, where 5 is a charged or neutral B meson, and 
the final state F involves mesonic states. Thus, one has to calculate {F\Hf.ff\B) , 
which involves the hadronic matrix elements {F\Oi{iJ,)\B). There are basically two 
competing approaches to calculating matrix elements: QCD factorization (QCDfac) 
and perturbative QCD (pQCD) 

The various WC's have been calculated at next-to-leading order fD^. It is found 
that C-i and Cs are negligible, so that, to a good approximation, the EWPs are purely 
(y — A) X (y — A). The consequences of this were investigated by Neubert and Rosner 
jS], and by Gronau, Pirjol and Yan (GPY) jB]. Using this approximation and flavour 
symmetries, GPY found two relationships between the EWP diagrams and the other 
{V - A) X {V - A) diagrams (T, C, E, A) in B ^ tttt and B ^ nK decays. 

From this brief summary, we see that there are many different aspects to hadronic B 
decays. In the present paper, we propose an approach which makes a connection to all 
of these. It is based on contractions. When calculating the amplitude for a particular 
decay, one must "sandwich" all operators of the effective Hamiltonian between initial 
and final states. All such terms have the form 

{qmqsqilqbbqeqrlqsb) . (5) 

(Dirac and colour structures are omitted for notational convenience.) Here gig2^3?4 
are the final-state quarks that eventually hadronize into the final-state mesons via the 
strong interactions. In our method, we are interested in the possible ways in which 
the final-state quarks can be produced in B decays through the effective Hamiltonian. 
These can be obtained by simply applying the basic rules of quantum field theory, and 
summing over all possible Wick contractions of all operators. 

Even though specific contractions between pairs of quarks in the decay process are 
involved, it is understood that any number of gluons may be exchanged between these 
quarks. Indeed, such exchanges are necessary in order for the final-state quarks to 
hadronize into mesons. The results of our general approach therefore hold to all orders 
in as- In our method we are not interested in the process of hadronization of the 
final-state quarks into mesons. If one wishes to calculate this, it is necessary to resort 
to models like QCDfac and pQCD. 

We will show that one can simply express all the diagrams in Ref. in terms 
of contractions. Since the set of all possible contractions follows directly from the 
effective Hamiltonian, this demonstrates explicitly the connection between diagrams 
and the OPE. As such, it constitutes a formal proof that the diagrammatic approach 
is rigourous. We also show that one reproduces the relations of GPY. Furthermore, 
we show that we can include gluons in the contractions order-by-order. (Of course, for 
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such an exercise to have any meaning, one has to argue that it is possible to calculate 
nonleptonic decays by including gluon exchanges order-by-order. The frameworks of 
QCDfac, pQCD and soft collinear effective theory (SCET) [7j have all discussed when 
and where an order-by-order expansion in as of hadronic decays is possible. ) 

Contractions were analyzed some years ago by Buras and Silvestrini (BS) [S]. At the 
most basic level, our method is simply a reorganization of BS. However, the connection 
between diagrams and contractions is much simpler in our approach. In addition, 
since the appearance of Ref . , we have obtained a variety of insights into nonleptonic 
decays. For example, there have been important developments in the various theoretical 
approaches. As such, the contractions can be separated into various classes, and the 
contractions in certain classes can be argued to be small because they arise only at 
0(l/mfc). In addition, numerous measurements in B decays have been made, and these 
help us to determine which contractions may be important. For instance, we know that 
certain rescattering contractions are in general smaller than other non-rescattering 
contractions. (We describe this in more detail in Sec. 2.1.) 

Using this general approach, we can express the amplitude for any decay fully in 
terms of contractions (indeed, this is what BS have done). However, as discussed 
above, some of these contractions can be argued to be negligible, either for theoretical 
or experimental reasons. We can therefore simplify the expressions with only minimal 
theory or model inputs. A fit to experimental data then allows us to obtain the various 
contractions. The fitted values can then be used to compare with theory calculations or 
make predictions for other processes. Note that our approach is valid not only for the 
SM but also for new-physics (NP) contributions. Hence, knowledge of the contractions 
from a fit to the data can give information about the nature of NP [9^ . 

In Sec. 2, we present the contractions method. The relation with the approach of 
BS is shown here. We make the connection to diagrams in Sec. 3. We show that one 
reproduces the GPY relations in Sec. 4, by working to all orders in a^. If one works 
only order- by-order in a^, one can make a connection between contractions and the 
matrix-element formalisms of QCDfac and pQCD. This is shown in Sec. 5. In Sec. 6, 
we examine whether there really is a "5 ttK puzzle." We conclude in Sec. 7. 

2 Contractions 

2.1 General Considerations 

We begin by describing the contractions in general. For a given B decay, one obtains 
many matrix elements of the form {qiq2q'i(li\(lbb qQq-j\q%b) [Eq. (0)], where Dirac and 
colour structures have been suppressed, is the B meson. This matrix element 
describes the decay B — > M1M2, where Mi^2 are mesons. The final-state mesons 
contain the quarks gi, g2, ^3 and ^4, but we are free to choose the quark assignments 
as we wish. In the following, we choose Mi = qiq2 and M2 = ^3^4. In addition, we will 
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write the final state for such B decays as 

1 



V2 



[M1M2 + M2M1 



(6) 



The reason we do this is as follows. In B —> nn decays, isospin requires that we 
symmetrize the final-state vr's. We choose the above form for the M1M2 final state in 
order to resemble B —>■ tttt. However, we stress that the analysis in this paper would 
not change with another choice, either of the Mi^2 quark assignments, or of the form 
of the M1M2 final state. For example, in 5 — > nK we can choose Mi = 51^2 = tt and 

= q^iQi = K, or Ml = qiq2 = K and M2 = q^Qi = tt without changing the physics. 
We refer to the freedom to exchange qiq2 ^ qsq^ as final-state symmetry. 

For a given b decay, there are 24 possible contractions: 



A = {qiq2q2,qM5bq&q7\qsb) 



r 



1 



1 



{qiq2qiqA\qbhq&q7\qsh) 



r 



1 



{qiq2q3q4\q5bq(iq7\q8b) 
I — I I — I I rn I 



G 

1= {qiq2q-iqA\qbbq%q7\q%b) 



r 



1 



M= {qiq2q-iqA\q5bqQq7\q?,b) 



r 



{qiq2qzqA\qbbqQq7\q8b) 
\ 



1 r 



mqsqAlqbbqeqrlqsb) 



r 



1 r 



{qiq2q3qA\q5bqeq7\q8b) 



B = {qiq2qm\q5bq6q7\q8b) , C 



{qmqsqilqbbqeqilqsb} , 



r 



1 



{qiq2qm\qbbq&q7\q8b) 



r 



.... "Fhn 
{qmqsqAlqbbqeqrlqsb) 

r 



H 

K = {qiq2q3qA\qrMf>q7\q?.b) 



1 



L 

L 



1 



{qiq2q3qA\qbbq&q7\q8b) , 



r 



1 



{qiq2qzqA\qbbq&q7\q8b) , 

I — I I I — I r-| — I 

{qmMAlqbbmrlqsb) , 



r 



1 r 



1 



r 



K= {qiq2q3qA\qbbqQq7\q8b) , = {qiq2qzqA\qbbq&q7\q8b) , 



Q 



r-hH n-hn 

{qmqsqAlqbbqeqrlqsb) 



r-F=Fr-F=hHn 



r= {qiq2q3qA\q5bq6q7\q8b) , U 



I — I I — I I r-| — I 

{qiq2q'iqA\q5bm7\q8b) , 



r 



1 



W. = {qiq2q3qA\qbbq&q7\q8b) , X = {qiq2qzqA\qbbq%q7\q8^) 



Not all the contractions are independent. For example, consider the two contrac- 
tions 



r 



1 



r 



1 



E{MiM2) = {qiq2q-iqA\qbbqQq7\q8b) , M{MiM2) = {q\q2q3qA\qbbq%q7\q8b) ■ (8) 

Here we label the contractions by (M1M2) with the understanding that Mi = qiq2 and 
M2 = q^qi- A contraction labeled by (M2M1) would assume the assignment M2 = qiq2 
and Ml = 53^4. Now, it is clear that the above contractions are not independent as we 
have M_{MiM2) — E_{M2Mi). This is just a consequence of the final-state symmetry. 
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Applying the same final-state symmetry to all contractions, we get the following 
equivalences between contractions: 

a b c d e f g h l i k l m n o p q r s t u v w x\ 
qroPkIktfJwvekcrabRkuLkg) 

Above, an element of the first line is specified by (M1M2) and is equivalent to the 
corresponding element in the second line labeled by (M2M1) and vice versa. Thus, 
any amplitude component can be be written in terms of the contractions A-X, or 
equivalently in terms of the "tilded" contractions of the second line. We therefore see 
that, of the 24 contractions, only 14 are independent. These correspond to the 14 
topologies of BS. We take the 14 independent contractions to be A, B_, C_, D_, E_, F_, G, 
IL, iL, K.1 L, N_, S_ and U_ (with or without tildes). 

Note that care must be taken in using these equivalences. Above it is assumed that 
the same operator is present in the two lines. If different operators are involved, these 
equivalences only apply if further symmetries are assumed (e.g. between Mi and M2). 
This will be important in Sec. 4. 

It is useful to separate the 14 independent contractions into four different classes. 
The advantage of this classification is that it becomes easy to implement theory inputs 
and certain classes of contractions can be argued to be small. 

• Class I: emission topologies: In these decays all the quarks in Hgff, apart 
from the b quark, are contracted with quarks in the final state. The spectator 
quark has no contraction with any quarks in i^e//- These contractions therefore 
have an "inactive" spectator quark. There are two such contractions, E_ and F_, 
which involve either an external or internal emission of final-state mesons. We 
therefore rename these with the ^EM^ label: 



I i I— t— 1 I I i I J j j 

E EM = {qiq2q-iqA\qbhq&q7\q^h) , F EMc = {qiq2q^qM^hqQq'j\qsh) . (9) 

The figures representing these contractions are given in Fig. ^ 

Note that in factorization we have EM ~ Fg^]^,f_^fj^,i^, where Fg_^j^j_^ is the B — > Mi 
form factor and J'm^ is the decay constant of M2. Hence, in this contraction M2 
is the emitted meson. Similarly, in factorization EM^ ~ Fg_j^j^fj^,j_^, and so 
here Mi is the emitted meson. This contraction therefore corresponds to colour- 
suppressed emission, and hence the subscript C is used in the label. 

• Class II: rescattering topologies: In these contractions the spectator quark 
does not contract with the quarks in Hgff and thus still remains "inactive." 
However there is a contraction between the quarks in H^ff. There are two possible 
contractions, involving the pair 55,^7 (contractions A and C_) or the pair ge??? 
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Figure 1: Emission topologies. Though not shown, the quark hues are dressed with 
gluons. 



(contractions G_ and K_). For each pair the quarks in the final state can have 
contractions between quark pairs belonging to different mesons {A and G) or to 
the same meson {C_ and K_). These latter contractions are expected to be small 
as they are Okubo-Zweig-Iizuka (OZI) suppressed. 

A possible A-tjpe contraction arises from the rescattering contribution of the 
tree operators. Since such rescatterings are usually referred to as penguin con- 
tributions, we rename this contraction P. Note that penguin operators can also 
produce this type of contractions. We rename the G contraction as PF since a 
Fierz transformation of the operator makes the G-type contraction look like an 
A-tjpe contraction. Finally, the OZI-suppressed contractions are renamed Pqzi 
and PFozi- 



] i— 1 I ] J I j ] i—i 1— 1 rn j 

A^ P = {qiq2q'iqimbqQq7\qsb) , G->PF= {qiq2q^qi\q5bq&qimb) , 

I — I I I I — I I — I I I — I rn I 

C_ Pozi = {qiq2q3qi\q5bq6q7\q8b) , K PFozi = {qiq2mA^^b%qj\%h) .(lo) 

The figures representing these contractions are given in Fig. |21 

In the QCDfac and pQCD approaches these rescattering contributions are per- 
turbatively calculable. For a given operator such contributions are suppressed 
by at least asijrii) relative to the contractions of the same operator belonging to 
Class I. For contractions involving charm quarks it is possible that rescattering 
may involve long-distance contributions that are not calculable perturbatively 
13 ^]. To get an estimate of the size of the rescattering contributions from 
charm intermediate states through h — »• ccs tree operators we can study charm- 
less B decays. Here there are no emission contractions for these operators; they 
can contribute only through rescattering contractions. Now, from experiments 
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Figure 2: Rescattering topologies. Though not shown, the quark hues are dressed with 
gluons. 

we know that these rescattering contractions are small and not of 0(1). If this 
were not the case, this would lead to too-large branching fractions for charm- 
less b s transitions such as B —>■ Kit, B —>■ (pK^ etc. The rescattering from 
b — i> CCS tree operators is then typically of the size of penguin amplitudes in 
charmless b ^ s transitions. Using flavour SU{3) symmetry, one can argue that 
rescattering through b ccd is also small. 

Hence, for a given decay and a given operator, the emission contractions of class 
I, which do not suffer any colour suppression, are generally larger than the rescat- 
tering contractions generated by the same operator. To take a specific final state 
as an example, we expect the colour-allowed decay B —> DsD to get its domi- 
nant contribution from emission topologies. Such arguments can also be applied 
to NP operators and have been used to argue for small NP strong phases from 
rescattering contractions jH]. 

As noted above, the OZI-suppressed contractions are expected to be small. In 
fact, the Pozi-type contraction can contribute to the decay B^ J/ipKs with 
a weak phase different from that of the dominant contribution through w-quark 
rescattering. If such contributions were significant, the result would affect the 
sin 2/5 measurement using this mode. The fact that the measurement JT] does 
appear to agree with predictions strengthens the claim that the OZI rescattering 
contributions are small. 
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Class III: annihilation/exchange topologies: We identify the annihilation 
and exchange topologies with the {B_, D_) and {H_, L_) contractions, respectively. D 
and L involve contractions between quarks in the same mesons, and are therefore 
OZI-suppressed. We therefore rename the annihilation contractions as A and 
Aqzi, and the exchange contractions as EX and EXqzi'- 



r 



B 



r 



1 



A = {qiq2q3q4\q5bq6q7\q8b) 
' ' 1 



r 



D ^ Aozi = {qiq2q3qi\q5bq6q7\q8b) 
These are shown in Fig. El 



H 
L 



EX 



I'?2g3'?4|g5&'?6'?7|g8&) 



EXozi{qiq2q3q4\q5bq6q7\q8b) 



(11) 




B-> A 




D-> A 

















H->EX 


b 


^ % 


1 1 
1 1 

1 77 1 









OZI 



L-> EX 



OZI 



Figure 3: Annihilation/exchange topologies. Though not shown, the quark lines are 
dressed with gluons. 

This class of contractions is suppressed by 0(l/m{,). Whether such contributions 
are sufficiently small is a matter of debate which will ultimately be settled by 
observing (or not) decays such as B^ K^K'^ or i?J] D^D^, which can only 
proceed through annihilation and exchange diagrams. 

The OZI-suppressed contractions in this class are again expected to be small. 
Their size can be obtained from the measurement of decays such as B^ — >■ J/ipcf) 
[T^ which can go through such OZI-suppressed contractions. 

• Class IV: topologies with annihilation/exchange + rescattering: This 
class includes contractions that have annihilation and exchange in combination 
with rescattering. These include iV, S_, J and £/. The N_ contraction can be 
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renamed PA for it is a "penguin-annihilation" contraction. Similarly, 5 is a 
"penguin-exchange" contraction: PE. Note that PA and PE are OZI-suppressed 
while the contractions U_ and J suffer an additional OZI suppression; we rename 
them as PAqzi and PEozi (Fig.|H): 



iv ^ Pv4 = {qiq2qzqMM(i(n\(isb) , s^pe = {qiq2qzq4\qM(iqi\qsb) , 

I — I I — I I Fn I I — I I — I I r-| — I 

J PAozi{qiq2qiqA\qbb%q7\%b) , H pEqzi = {qmrn^l^bmilqab) , (12) 




Figure 4: Topologies with annihilation/exchange + rescattering. Though not shown, 
the quark lines are dressed with gluons. 

The contractions in this class are expected to be tiny as they are both 0(l/m5)- 
suppressed and OZI suppressed. 

To summarize: there are 14 independent contractions. These can be separated into 
four classes: {EM, EMc}, {P, PF, Pozi, PFozi}, {A, EX, Aozi, EXozi}, 
{PA, PE, PAozi, PEqzi}- In what follows, we will almost always express amplitudes 
in terms of these contractions. However, we will occasionally write amplitudes in terms 
of the 24 contractions A-2L [Eq. 0]. When we do this, we will be sure to give a warning 
of this fact. Thus, unless there is an explicit statement to the contrary, the reader can 
assume that all amplitudes are given in terms of the contractions EM, EMq, etc. 

Buras-Silvestrini (BS) also have 14 topologies, described by rather complicated, 
nonstandard figures. They label these as DE (Disconnected Emission), CE (Connected 
Emission), DA (Disconnected Annihilation), CA (Connected Annihilation), DEA (Dis- 
connected Emission- Annihilation) , CEA (Connected Emission- Annihilation) , DP (Dis- 
connected Penguin), CP (Connected Penguin), DPE (Disconnected Penguin- Emission) , 
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CPE (Connected Penguin- Emission), DPA (Disconnected Penguin- Annihilation), CPA 
(Connected Penguin- Annihilation), DPA (Disconnected Double- Penguin- Annihilation), 
CPA (Connected Double- Penguin- Annihilation). Their eight OZI-suppressed topolo- 
gies are DEA, CEA, DPE, CPE, DPA, CPA, DPA, CPA. As such, the starting point 
of our method is very similar to theirs. However, what they do is rather different from 
what we do in this paper. 

BS's main aim is to produce a manifestly scheme- and scale- independent formalism 
for hadronic B decays. In the subsequent sections we show that this method reproduces 
results which have been obtained using other techniques. In particular, the results of 
diagrams, GPY, etc. can be simply expressed in terms of contractions. 

BS stress the need to include the "small" OZI-suppressed contractions in order 
to obtain scheme- and scale-independent results. We agree with this observation, in 
principle. However, the BS paper was written before the observation that certain con- 
tractions are suppressed by factors of 0(l/mfe). Thus, in practice, these contributions 
are indeed small, and it is reasonable to neglect them when computing the amplitudes 
for hadronic B decays. Including the small contractions makes the expressions for the 
amplitudes unnecessarily complicated. 

Finally, BS show how to write diagrams in terms of their contractions (Sec. 10 of 
their paper). However, these diagrams are not those described in Ref. and used 
elsewhere. In the present paper we show how the diagrams of Ref. can be written 
simply in terms of contractions. 



2.2 B ^ttK Decays 



In this subsection, we show how to compute the possible contractions for a given 
B — > M1M2 decay. We do so by considering the four B — > ttK decays: B^ n^K^, 
B- ^ n^K-, 50 ^ ir+K-, 5° ^ n'^K''. 

We begin with B~ —>■ tt~K^. The operator Tix that generates tree amplitudes 
involves Oi and 02- Then 



7^ E ((--^° 



B-) + (K^-R- 



B- 



contractions 



V2 



-0\c 



12 ' 



where 

For 01: 
(T- 



contractions 



B- 



TT 



B- 



0\u 
)l 



71 E {{^dds 

V ^ cont'i 



uhsu \ub) + Idsud uhsu \uh) 



contractions 
lu 



)= [(A-(7r-ir°) + 5r(7r-ir°)) + &';{K\-)+U:{K\~)) 



V2 



(13) 
(14) 



(15) 
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In the above, we have expressed in the amphtude in terms of the contractions A-X of 
Eq. ((7j). Here and below, the primes on the contractions indicate a 6 — s transition. 

Now, final-state symmetry indicates that (^°7r~) = A^^{7i~K'^) and^;^ {K^7r~) = 
B^i{7r~K^). Using this, and changing the notation to that used in the four contraction 
classes, we then obtain for -k^K^ (dropping the (M1M2) label): 

{T-y^ = y2(P;" + A'^) . (16) 

(Here, P' and A' indicate penguin and annihilation contractions, respectively.) In 
the above we have shown how final-state symmetry is used to include the Q and R 
contractions. From here on, final-state symmetry will be assumed and not shown 
explicitly. 
For 01: 



(T °)J = — p ^ [{udds cbsc\ub) + (dsud chsc\uh)^ . (17) 

V ^ contractions 

Only the P' contraction is allowed and we get 

(T-°)5 = V2Pf . (18) 
For O2 and O2, only the colour structure is different, so that 

(T-O)^ = v^(P- + A-), 

{T-% = V2P^\ (19) 



givmg 

T-° = V2ci [Xi'^ (/^'" + A'n + Xi'^P^'] , (20) 

where a sum over z = 1, 2 is understood. 

The tree pieces of the other three B ttK decays can be found in the same way. 
In all, we have 

T-° = V2c.[xi^\p:-+A'n+xi'^pn, 

T°- = Q l-Xl'^ [P'^ + A'l' + em;" + EM'^;) - Xi'^P^'] , 
r+- = V2Q[-A(f)(i^'" + EMf) - Xi'^P^'] , 

T0° = c,[Xi^\Pt-EM'^^) + Xi^'>Pn. (21) 

Note that the quadrilateral isospin relation for B ttK decays is respected: 

T-° + V2T^- = T+- + V2T°° . (22) 

We now turn to the EWP operators. We make the approximation of neglecting C7 
and Cg, so that the EWP operators are purely {V — A) x {V — A) and are given by O9 
and Oio alone. 
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In these operators, one sums over the hght quarks, multiphed by their charge. Below 
we consider g' = u,d,s (with e„ = 2/3, Cd = Cs = —1/3). Below we will always assume 
isospin symmetry, so that the contractions of u and d quarks are equal. However, we 
do not assume flavour SU(3) symmetry here (under SU(3) the contractions of u, d and 
s quarks are equal). Instead the u/d and s contractions are labeled individually. 

The EWP contractions for the four B — > nK decays can be found in the same way 
as was done for the tree operators. We have 



p-O 

EW 




-ipF: - '^EX', + ^Em'^, + ^ {Pt + PF'i 


) , 




pO- 

EW 


- ^A^'V 

- ^\ Q 


^PF[ + EM[ + hx[ + hn'^, - i (P:^ + PF'n 




p+- 

EW 




\PF[ - \eX', + \em',, - i {P[^ + PFn 






pOO 

EW 


- ^\ Q 


-^PF[ + EM[ + ^EX[ + ^EM',^ + ^iPt + PF'n] ■ (23) 



Here, a sum over i = 9, 10 is understood. Note the appearance of the terms EM- and 

EM^^. We remind the reader that the tilde indicates that the final state is M2M1 
as opposed to M1M2. Thus, for B — > ttK decays, the EWP contractions involve final 
states of the form Ktx. (To be more precise, EWP contractions involve ttK terms which 
are equivalent to Kt: contractions, as per the equivalence table following Eq. (jS)).) 

In the above, w-quark and rf-quark contractions are taken to be equal (isospin 
symmetry). The u/d contractions have no label, while the contractions of s-quarks are 
labeled by an 's' superscript. Because the operators in the two sets of contractions are 
different, the quarks are in different positions, and the equivalences due to final-state 
symmetry do not apply here. On the other hand, in the SU(3) limit, w-quark, (i-quark 
and s-quark contractions are equal, and the effect of SU(3) symmetry will be explicitly 
worked out in Sec. 4. One sees that, once again, the isospin quadrilateral is respected 
above. 

We note in passing that there are rescattering contractions involving the c quark for 
the EWP penguin operators which are represented by -Pg'^io ^^"^ -^-^g'lo- However these 
contributions are expected to be tiny relative to the rescattering from the tree operators. 
This is unlike the case of rescattering contractions involving the u quark where the EWP 
contractions are enhanced by CKM factors relative to the tree contractions, so that 
both contractions are of the same size. 

Finally, we turn to the gluonic-penguin operators. Here no WC's can be neglected, 
so that all operators 0^-0^ must be included (i.e. we take into account both {V — A)^ 
{V — A) and iy — A) x [V + A) Dirac structures). We have 



p-O 
pO- 



-V2\ 



2PF', + PFf + EX', + EM^i + P- 



Ci 



2PF', + PF'f + EX', + EM^, + P'^ 
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pOO 



V2\rc, 2PF[ + PF'l + EX[ + EM^, + P^' 



V2\['^^^2ci hPF'. + PF'f + EX' + EM' , + P, 



(24) 



which respects the isospin quadrilateral relation. The sum is over i = 3, 4, 5, 6, and the 
index indicates the Dirac structure. 

In the above, we have made no assumptions about colours. That is, the quarks can 
have any colour, representing the exchange of any number of gluons. Thus, the above 
equations hold to all orders in as- (We will work order-by-order in Sec. 5.) 

In this subsection, we have shown how to write the matrix elements for the four 
B —>■ ttK decays in terms of contractions. This can be done for any hadronic B decay. 
(Indeed, this is essentially what BS have done.) One important point is that, while we 
have expressed matrix elements in terms of contractions, we have not evaluated these 
matrix elements. This requires an additional method, such as QCDfac or pQCD. 



3 Connection to Diagrams 

All B ttK decays receive contributions from the penguin diagram, P'. (As with 
contractions, a prime on a diagram indicates a.b ^ s transition.) This diagram actually 
contains three pieces, corresponding to the identity of the internal quark: 

^' = Ai^)p: + AWp^ + Af)p;. (25) 

Below we use this form for P'. Note that the unitarity of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix has not been used. This is because the unitarity of the CKM 
matrix was not used in writing the effective Hamiltonian [Eq. (jT))]. 

When writing the B ttK amplitudes in terms of diagrams, it is conventional to 
absorb all A^''^ {p = u, c, t) factors into the diagrams themselves. We then have 

= [pL + aU[Pc] + [pI-IpLI- + IpLI + Ip'/w], 
v/2A°- = [~r -c -p^-A'] + [-p;] 

2 2 1 

_L \_p' _ p' _ _p'C _ Zp'E _ tp'Pn] (Of.) 
' [ t ^ EW g-* EW g-* EW g EW i 5 l^^J 

A^- = [-r-PL] + [-p:] + [-p;-lpz^ + lp'z^-lp'i^], 

v^^°° = [-c' + pL] + [p:] + [p;-pL-~IpL'w-Ipe''w + Ipe''w]- 

In Ref. in] it was shown that there is an EWP diagram corresponding to each of the 
T, C, Pu, A, E, PAu diagrams. These are included above. The first, second and third 
terms in square brackets are proportional to A^'^^ A^'^^ and A^'''', respectively. 
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In the previous section it was shown how to write the contributions from tree, EWP 
and gluonic-penguin operators in terms of contractions. Using the above expressions, 
we can now write diagrams in terms of contractions. The tree-operator contribution 
is given in Eq. (pT|). Comparing with the first and second terms in parentheses above, 
we obtain 

r = V2Xi'^CiEMl'' , 

A = v^Ai^^QAf, (27) 

where the sum over i = 1,2 is still understood. This matching looks very natural. 
Graphically, an EM' emission contraction really looks like a T' diagram, and similarly 
for the colour-suppressed emission contraction, EM'^, and C (Fig. Q). and P^ are 
both described by the same type of contraction (Fig. |21), which represents rescattering 
from the operators Oi^2- Finally, the A' contraction can be identified with annihilation 
diagram (Fig. El). 

The contribution from EWP operators is given in Eq. Comparing to the EWP 
terms above, we have 



p' = 

EW 




pic _ 

^ EW 




pIPu _ 

EW 


-v^^Ai^)Q[(ppf-pp:)+/^- 


plE _ 
^ EW 


-V2h['^c,EX[. 



Finally, the gluonic-penguin contribution is given in Eq. 



This implies 



P; = -V2\Tc, 2PF' + PF'f + EX' + EM,, + P' 



{2i 



(29) 



Now, we note that Eqs. ()27|) . ()28|) and (j^^ above are missing some diagrams: E, 
PAu^ct, Pew s^nd Pew""- Order to obtain these, we must use a different decay. Here 



we choose the three B ^ tx-k decays: B' 



B" 



TT ' TT 



B^ 



7r°7r° (recall 



that isospin symmetry has been assumed, so that the placement of the final-state pions 
is unimportant). 

We simply present the results here; they can be derived using the techniques de- 
scribed previously. The tree contractions are 

T"""° = -Af q(EM;' + EM^J , 

j^+n- ^ ^^{d) ^^^u + p." + EX^ + PE^) + A^'') (P^^ + PE^)] , 

T^°-" = c, [Af (-PM^, + i^" + EX^ + PE^) + A^'^) (i^^ + PE^)] . (30) 
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The EWP contractions are 
3 



T^TT TT 



\f^c,{EMi + EM, 



EW 



CiJ : 

'^EMci + ^ (A + PFi + PAi) 

-l{P^ + EX, + PEi) - i {PFt + PAD 
EMi + ^EMai - ^ (A + PFi + PAi) 

+ l{Pi + EXi + PEi) + i {PFI + PAf) 



(31) 



The contractions from the gluonic-penguin operators are 

rrO 







= V2Xi:'^Ci \(2PF, + PFI + EX, + EM a, + P, 



— '^t 



+ {2A + 2PA + PAt + PEi)] 
(2PFi + PFI + EX.; + EMci + P) 

+ (2A + 2PA + PA' + ^^0] 



(32) 



In the above, the absence of a prime on the contractions indicates a b ^ d transition. 
Note that, in all cases, the triangle isospin relation for S — > tttt decays is respected: 



(33) 



In terms of diagrams, the amplitudes for the three B ^ tttt decays are given by 



^7r+7r~ 



-^[(T + C) + (P,,. + P,^J] , 

-[(T + P^ + E + PA,,) + (Pe + PA,) 



-l-(P-l-PA-l--P'^ -1-1 _1p^ _1 _ 1 pPAu \ 

"r V t t 3 3 3 3 3 



[{C-Pu-E- PA„) + (-P, - PA,) 
+ (-P, - PAt + P,^ + ^P^^ - ^P^^ + ^Pf 



(34) 



1 P^" _L 1 pPAu\ 
EW ' 2 2 / 



Comparing the amplitudes in terms of contractions and diagrams, we see that the 
missing diagrams are given by 

E = y/2\^^'^CiEXf , 
PA„ = y/2\^^CiPE^ , 
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Table 1: Connection between diagrams and contractions. For the tree diagrams (T, C, 
Pu, Pc, A, E, PAu, PAc), a sum over i = 1,2 is understood; for the EWP diagrams 
(Pew, Pew, ^ew, Pew, Pew, Pew), the sum is over i = 9, 10; for the gluonic-penguin 
diagrams {Pt, PAt), i = 3, 4, 5, 6. q is d or s, depending on the b decay. 



Diagram 


Contraction 


Diagram 


Contraction 


T 


V2Xi^kiEM^ 




p 

EW 


-V^IXfciEMi 




C 


V2\^^h,EM-, 




pc 

EW 


-V2lXi'^c,EMa 


i 


Pu 






pPu 
^ EW 


-V2-IX\'^c,[{PFt-PF,) + Pi] 


Pc 


-V2\^^h,P^ 










A 






pA 
^ EW 


-V2fA(^)c,A, 




E 






pE 
^ EW 


-V2lXi''^CiEXi 




PAu 






pPA-u 
EW 


-V2lx't''''c,[{PAt-PA, 


) + PE^\ 


PA, 


-V^X^^hiPE', 








Diagram 


Contraction 




Pt 




Ci 


2PF, + PFt + EX, + EMci + P! 




PAt 


-V2X 


^Si [^Ai + 2PA + PAl + PEi] 






PA, 


= V2X^f^ 


cPEt , 








PAt 


= V2X 


id) 


Q {2A, + 2PA, + PAl + PEi) , 






pA 

EW 


2 


xi ^CiAi 







PZ- = - V2 ^A^'^^Q [{PAl - PA,) + PEt\ . (35) 

For diagrams which appear in both B — ^ t:K and B ^ Tin decays, the expressions for 
the connection to contractions are the same as in Eqs. (jTfj) . and (P^ . except that 
the prime is removed. 

The connection between diagrams and contractions is summarized in Table ^ One 
sees that the expressions for Pt and PAt are quite different from those for P^ c and 
PAu^c, respectively. This is because Pt and PAt are derived from contractions involving 
O^-Oq, while P„ c and PAu^, involve Oi and O2. 

Note that this connection involves only 8 of the 14 independent contractions: EM, 
EMc, PF, P, A, EX, PA, PE. The six OZI-suppressed contractions Pozi, PFqzi, 
Aqzi, PAqzi, EXqzi, PEqzi are not used. The reason is the following. Each of these 
six involves the contraction of qi and q2 and/or gs and ^4. This can only happen if the 
final state contains neutral particles of the form qq. In the decays we have considered, 
only 7r° obeys this criterion. However, since 7r° = [dd — uu) / \/2, and since we have 
assumed isospin symmetry, the six contractions vanish. It is only if we consider final 
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states involving a <^ (ss) or J/\l/ (cc) that these contractions will contribute. 

4 Connection to GPY 

In this section, we show that, using contractions, one can reproduce the main results 
of Ref. in], by Gronau, Pirjol and Yan (GPY). GPY use a group-theoretical formalism 
to analyze B ttK decays. On the other hand, the analysis using contractions is 
more direct. Throughout Ref. GPY neglect the Wilson coefficients (WC's) C7 and 
Cg, so that the EWP operators are purely {V — A) x {V — A). In what follows, we 
will make the same approximation. GPY also assume flavour SU(3) symmetry. We 
will (eventually) do likewise, but we will show how this assumption is necessary in our 
approach. 

The first result is the following. The WC's obey Ci/c2 = cg/cio to about 5%. In 
the limit in which this equality is exact, GPY note that the EWP amplitudes (-P^w 
P'f^, P'^'l, P'f^, P'/^-) are proportional to the tree operators (T', C, P^, A', E', 

To see how these relations emerge in our formalism, we first concentrate on the 
diagrams T' and P'j^^ in i? — ^ nK decays. From Table [TJ we see that these are propor- 
tional to EM and EM, respectively. These imply that the final states ttK and Ktt, 
respectively, are produced. Now, we note that for P^^^,, only the up-quark piece of the 
electroweak operators contribute - this is denoted with an index u. That is, P^,^ in- 
volves only Oq^q ~ sbuu (as usual, the colour and Dirac indices have been suppressed). 
However, in the SU(3) limit, s-quark contractions are equal to u-quark contractions. 

Thus, ~ {Saba)v-AiU(3Uf3)v-A = {Uaba)v-A {s pUf3)v-A Equivaleutly, IT = K, 

so that the final states ttK and Ktc are the same. Thus, EMg = EM[. Things are 
similar for Oio and O21 so that EM^q = EM^. We therefore see that 

r = V2\l'\ciEM[ + C2EM!,) = V2\l'^ci (eM[ + '^EM'^ , (36) 

while 

P',^ = -V2h['\c,m', + c,,EM[,)'''S'^ -V2hfc{EM[ + ''f^EM'^ . (37) 

In the limit where Ci/c2 = cg/cio, P'^yy is proportional to T'. 

The argument is much the same for C and P^'^. In the SU(3) limit, we have 

EM'c9 = EM'c^ and ^'^lo = EM'^^. Thus 

C = V2Xi'\c,EM'^, + C2EM'^,) = v^A(^)ci (^EM'^, + "fEM'^^ , (38) 

and 

p'-^ = -V2h['\c,m',,+c,,m',,,)''i^'^ ^ . (39) 
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As above, is proportional to C in the limit where Ci/c2 = Cg/cio- 

Things are simpler for the other four EWP diagrams P'/^, P'ewi P'ewi P'^w" ■ 
In the SU(3) hmit, PFf - PF'^ = 0, so that P'/^ is just proportional to P^ {i = 9, 10). 
Similarly, PA'^'-PA[ = in the SU(3) hmit, and PJ^" is proportional to PE[. Thus, aU 
four EWP diagrams involve the same contractions as the corresponding tree diagrams 
P^, yl', E' and PA'^- For P^ and P^^ we can therefore write 

P'u = V2A W {c,P[ + c^P^) = y2AWci (p; + ^P^) , (40) 

and 

P'l^ = -^^A(^^(c9P; + CioP^) = -V2^aJ^)c9 (p; + ^P,') . (41) 

In the limit where Ci/c2 = Cg/cio, P^^ is proportional to P^. The argument is identical 
for the other three EWP diagrams P^^, P'^w^ Pew"" diagrams A', P', PA'^. 

We therefore see that the proportionality constant between the EWP and corre- 
sponding tree diagrams is —{?>/2){cg/c\)\^''/\^\ which agrees with GPY. In particular, 
we have 

P'ew 3 \ ^ P'ew 3 Xt Cg 



(42) 



The second result is two relations between EWP and tree diagrams: 

P,^(p-^rf°) + V2P,^(p-^7r°i^-) = ^^l±^R(C' + r) , 

2 C1 + C2 

P^^iB'' ^ TT+X-) + Pw(S- ^ 7r-X°) = 

3 c9-cio ^^^, + C" - T' - P) - |^^±^P(A' -C'-T' + E) , (43) 

where R = aI'V^I"^- All the dia grams in the above relations are for B — > tiK decays, 
with the exception of P, which is for B — > tttt. However, since flavour SU(3) has been 
assumed, B 71K and B ^ rni diagrams are equal. 

In order to reproduce these relations, we make the following observation. In the 
limit of neglecting C7 and Cg, both EWP and tree operators are {V — A) x {V — A). In 
this case, it is possible to do a Fierz transformation to exchange the position of and 
Qq without changing the Dirac structure of the operators. This results in Oi ^ Oio 
and 02^09. If wc also take the SU(3) limit, in which case one can switch u and s 
quarks, wc have Oi = Og and O2 = Oiq. Given the equivalence of different operators, 
this implies that certain contractions are pairwise the same (within each of the four 
contraction classes): 

PMi,9 = EMc2,io , EM2,io = EMci,9 , 

Pl,9 = -P-^2,10 ) P2,10 = PFl,9 , 

-foz/1,9 = PFozi2,10 -I Pozi2,10 = PPoz/1,9 , 
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^1,9 = -£'-^2,10 , ^2,10 = £■-'^1,9 , 

= EXozi2,10 , Aozi2,W = EXozil,9 , 

PAi^Q = PE2,io , PA2,io = PEi^Q , 

PAozil,^ = PEozi2,10 , -P^oz/2,10 = PEozil,9 ■ (44) 

Consider now the first relation. The left-hand side of the equality is 

i=9,10 ^ 

= V2^Xl'^ {cgEMg + CioEM[q + CqEM'^^ + c^^EM'cio) 
^^'^ ^/2-\['^ (cgEM[ + c^oEM^ + CgEM'^, + c^oEM'^^) 
Fierz+SU(3) ^\is) ^^EM[ + c^^EM'^ + c^EM'^ + C^oEM[) 

= ^\\^\c^ + Cio)(SM( + EM'^) . (45) 

The right-hand side is 

i=t2 2 Ci + C2 

= ^3^W C9 + Cio ^^^^^; ^ ^ CiEM;,! + C2^M^2) 

2 Ci + C2 

F.e..+5f7(3) ^3^W^9 + C10 ^ ^ ^ 

2 Ci + C2 

= V2^A«^^±^[(ci + C2)£;M( + (ci + C2)£;M^] 

2 Ci C2 

= \/2^ aJ^) (c9 + Clo) (^M( + EM^) . (46) 

We therefore reproduce the first GPY EWP-trcc relation using contractions. 

We now turn to the second relation of GPY. The left-hand side of the equahty is 

^ Ci^"^ [em'^, - EXl) = x/2|aJ^) (c9^c9 + cio^cio - Cg^^g - CioEX[,) 



1=9,10 



= V2^A;^^ {cgEM^ + cioEM[ - cgEX[ - c^EX'^) (47) 



and the right-hand side of the equality is 



-3 



=1,2 4 



Y: V2^^X^^ Q^^^(^: + EM',, - EM[ - EX'^ 

^ ^ L Ci - C2 
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^C2 



- - EM',, - EM[ + EX[] 

Ci + C2 
4 Ci — C2 

-Ci^M( - C2^M^ - ciEX[ - C2EX'^) 

4 Ci + C2 

-CiEM( - C2EM'^ + ciEX( + C2EX'^) 

^3^(,) C9-Cio _^ _^ ^ 

4 Ci - C2 

-ciSM{ - C2^M2 - ciEX[ - ca^X^) 
_y23^(^)^^±^ (c^^X' + C2EX[ - ciEM' - c^EM'. 

4 Ci + C2 

-ciEM( - C2EM'^ + ciEX( + C2^X^) 

^3^(,) C9-Cio _ _ ^ _ 

4 Ci - C2 



v^^a1^^^^^(ci + C2)(i5;X^ + EX[ - EM', - EM',] 

4 Ci + C2 



= 72^ aJ^^ (cgEM^ + cio^M{ - C9EX( - cio^X^) . (48) 

This proves the second EWP-tree GPY relation in — ttK. 

We have therefore shown exphcitly that, using contractions, we reproduce the main 
results of GPY jH] for i? — > nK decays. 

It is worth making one final remark here. We have mentioned in Sec. 2.1 that the 
contractions A'^ and EX'^^ are higher order in (l/rrife). If these contributions are indeed 
small, then the second GPY relation can be separated into two relations, one leading 
order (LO), the other next-to-leading order (NLO). These contain the large and small 
contributions, respectively. We have 



LO 



3 cg-cio^^^, - T') - -^-^±^R{C' - T') , 

4 Ci - C2 ^ ^ 4 Ci + C2 ^ ^ 

J JVLO 

4 Ci - C2 ^ ^ 4 Ci + C2 ^ ^ ^ ^ 

Of course, each of these relations is individually satisfied. 
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5 Connection to QCDfac and pQCD 



All the results in the previous sections hold to all orders in ag- As mentioned earlier, 
in our general approach, we are interested only in the possible ways in which the final- 
state quarks, which evolve into the final-state mesons, can be produced in B decays 
through the efi^ective Hamiltonian. That is, we do not attempt to calculate the strong 
interactions of the final-state quarks through the exchange of gluons and their ultimate 
hadronization into the final-state mesons. However, it is also possible to work order- 
by-order in Us by including the gluons explicitly. This is useful as it allows us to make 
a connection with QCDfac and pQCD, both of which work ordcr-by-ordcr. They do 
this because in their frameworks the nonlcptonic amplitudes are expanded in powers of 
as{nh) {iJ-h ~ \J-^QCD'tTib) for gluon exchanges involving the spectator quark, or as{mb) 
for the remaining gluon exchanges. 

The first calculation involves gluon exchanges only between quarks that belong to 
the same meson. In our approach, such gluon-exchangc effects arc taken into account 
automatically, and so we call this the 0{a^) ("zero-order") calculation. This is in fact 
just the factorization assumption in calculating nonleptonic amplitudes. 

The exchanged gluon between the quarks in a meson can be soft and hence these 
contributions are not pcrturbatively calculable. In the pQCD approach, the soft gluon- 
exchange contributions involving the soft spectator quark in the B meson are argued 
to be highly suppressed. Thus, the soft spectator quark has to be boosted up in 
energy through gluon exchange in order to end up in the final-state light meson with 
energy E ~ M^- Since the ghion is relatively hard such contributions are calculable 
pcrturbatively in pQCD. However, it is not entirely clear that soft-gluon exchanges 
involving the spectator quark can be neglected, and in the QCDfac approach the soft- 
gluon contributions are absorbed in physical form factors. 

In the second step we consider gluons exchanged between quarks belonging to dif- 
ferent mesons [0(a])]. In our approach we do not calculate the amplitudes explicitly. 
Instead, we work out the general structure of the contractions when one- gluon ex- 
changes are included among the quarks that emerge from the B decay through the 
weak effective Hamiltonian. 

In the pQCD approach it is argued that such corrections are calculable pcrtur- 
batively, while in the QCDfac approach the gluon exchanges involving the spectator 
quark are sensitive to long-distance physics. Thus, their effects, like the form factors, 
are represented by additional unknown hadronic quantities that may be obtained from 
a fit to data. 

5.1 Order zero 

In this section we will consider the zero-order calculation. Prom here on, we concentrate 
only on T', C, P^,^, and P^*^ in P nK decays. In this case, with no gluons between 
quarks of different mesons, it is possible to take colours into account simply by counting 
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them. For example, for T we have 



T 



i=l,2 



r 



1 



r 



1 



= Ci{qi^q2xq'iyqiy\qbJ>aq(Sfiq7p\q%zhz) + C2{qixq2xq3yqiy\q5ab(3q6(iq7a\q8zbz) -(50) 

In the above, the subscripts a, /5, x, y and z are colour indices. Since gluon exchanges 
among quarks belonging to different mesons are neglected, the contractions are simply 
proportional to delta functions in the colours: 



r 



1 



r 



1 



i=l,2 



CiEMi = Ci6xz5xaSyfsSyfs6az{qiq2q3q4\q5bq6q7\q8b) + C2SxzSxa5ya5yf36f3z{qiq2q3q^\q5bq6q7\q8b) 



r 



1 



r 



1 



CiN^ {qiq2q?.qi\qM&q7\q8b) + C2Nc{qiq2q3qi\qM&q7\q8b) 
CiNl^M + C2N^'EM 



C2 



(51) 



where the bar on EM is added to stress the fact that colour effects are extracted. The 
idea is the following: a priori the matrix elements EMi and EM2 are different. But 
since the difference between Oi and O2 is only their colour structure, once we have 
taken this into account by explicitly adding the effect of colour, the remaining matrix 
elements are identical. Doing this for T', C, 
following relations 



and P'^yy, it is easy to derive the 



r 
c 



p' 

EW 



v^(ci + ^)Nl EM' 



V2{^ + C2)NlEMa' 



V2{cg + ^)N'^EM 



N, 



^V2{^ + c,o)N^EM',. 



(52) 



With SU(3), and taking into account Fierz transformations, we have EM' = EM'^ = 
/ / 

(which hold to 



EM = EM^. Then we obtain the two EWP-tree results of Eq 
all-orders), and a third relation: 



V 



(ci + CsA^c) 
(ciiV, + C2) 



0.17±0.1 



(53) 
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The numerical value of this ratio is obtained as follows. At leading order the WC's 
take the values ci = 1.081 and C2 = —0.190 (these values are taken from the last paper 
in Ref. [S]). Thus, taking Nc = 3, there is some cancellation between the factors in the 
numerator of Eq. (j^Hj) . Now, the WC's are calculated using a renormalization point of 
/i = rrib- However, in fact we do not know the precise value of /i - all we know is that 
it is 0{mh)- Allowing the values of the WC's to vary leads to an uncertainty on the 
above ratio. The central value uses the above values for ci and C2, but the (estimated) 
error corresponds to allowing /i to vary. The bottom line is that the precise value of the 
ratio C /T' is uncertain. Still, although this ratio can be zero in the case of complete 
cancellation, there is an upper limit of approximately 1/3. We note in passing that the 
naive estimates of colour suppression are supported by the WC calculations. 

Above, we have concentrated on T', C", P^^y and P'^yy in 5 — >• tiK decays. However, 
one can perform the above procedure for other contractions/diagrams. If one does so 
for the annihilation and exchange diagrams A and E, one finds the following interesting 
relation. At O(a^), the matrix elements are equal, so that 



A 



Ci + ^ 



6.0 . (54) 



We therefore see that, at leading order, the diagram A' is expected to be much larger 
than E'. BS found a similar result, but in the large Nc limit. 



5.2 Corrections with 1 gluon 

We now turn to the 0{al) QCD corrections. We introduce a single gluon between all 
pairs of quarks in each diagram/contraction. We illustrate this procedure by consid- 
ering the ii^M-type contraction for tree operators in B ^ ttK (Fig. E)). There are ten 
possibilities for the gluon, shown in Fig. El Note that these include a possible gluon 
exchange between the quarks within the same meson - we will be able to absorb these 
corrections into the order-zero result. 

We introduce the following notation: Zij is the j^^ QCD correction (j = 1, 10) 
to the Z-type contraction {Z = EM, EMa,..., ) of O,. Zjj is the same, except that 
the colours have been explicitly extracted: Zij = {colour factor) x Zij). 

The idea is the same as for the 0{a^) effects: the only difference between Oi 
and O2 is the colour structure. For example, EMi^i and EM2,i differ only in their 
colour structure, and so their matrix elements are equal once the colour effects have 
been extracted. The problem is more complicated here because there are several QCD 
corrections, e.g. we cannot relate EMi^i to EMi^2- Thus, the overall matrix elements 
are not equal unless all contributions appear in the same linear combination. We will 
see this explicitly below. Note also that matrix elements are still not evaluated. 

We begin by computing the colour factor for EMi^i 
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Figure 5: EM' contractions of tree operators for B ttK. 




Figure 6: One-gluon corrections to EM 



+ ^ = . 



2Nc 2 

Then EMi^i vanishes by colour arguments. For EM2^i, 

^yx^xm^ /3y^ zfB^nz 



1 



2N^ 2 2 I m 



The total contribution of the correction 1 to the contraction EM is then 



K 



(55) 



(56) 



(57) 
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A similar calculation has been performed for the remaining nine QCD corrections 
of the contraction EM , with the following result: 



io / 1 \ 6 



+ Ci + 



Nr. 



1 



^^c / i=7 



10 



(58) 

Because the coefficients of the two terms are not the same, the overall 0{al) matrix 
element of Oi is not the same as that of 02- 

We can apply this to the T' diagram in i? ^ -kK decays. Let T' = T'° + T'^ where 
T'" is the piece without gluons [0(a;°)] and T'^ is the piece with one-gluon corrections 
[0{al)]. We have 



T 



•II 



V2 



C2- 



1 - 



^^c / i=l 



Cl + 



C2_ 

Nr. 



N^ 



10 



1 - 



.(59) 



However, an examination of Fig. IHl reveals that the corrections i = 7,8,9, 10 all cor- 
respond to gluon exchange between quarks within the same meson. As such, these 
corrections can be absorbed in the 0{a^) result. We will therefore rewrite Eq. EHlas as 



rp>0 
rpll 



v^(ci + ^)NlEM', , EM',^EM' + ^[l- ^EM, 



(60) 



This procedure can also be carried out for the EMq contraction. Fig. [7| shows the 
E'Mc-type contraction for colour-suppressed tree operators in i? — ^ ttK. As before, we 
consider a single gluon exchange between all pairs of quarks. There are ten placements 
for these gluons, shown in Fig. |H1 

Taking all one-gluon QCD corrections into account, we find 

10 N"^ { l \ ^ / ci \ N^ ( 1 \ ^° 

Y.{c,EMa,,^+C2EMa,..) = ci^ ( 1 - ^1 E^^c.+ (^ + "V 1^ ( ^ " ]^ ) E™c. 

(61) 

Thus, as for the EM contraction, the overall 0{a\) matrix element of Oi is not the 
same as that of 02- 

We now apply this to the C diagram in i? ^ tiK decays: as with T' we define 
C = C'o + C'l, with 



V2{ciNr + C2NI) EM[ 



V2 



N? 



C ' 
6 



^^c / 1=1 



-i^(1-1^1E™o.+ (^ + C2)^(1-^)E™ 



Nr. 



N^ 



10 



N^ , 

^^c J i=7 



(62) 
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Figure 7: EM'^ contractions of colour-suppressed tree operators for B tiK. 
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Figure 8: One-gluon corrections to EMc 



K 



However, like the T' contribution, the corrections i = 7, 8, 9, 10 can all be absorbed in 
the 0(a°) result. We therefore rewrite 



C 



(63) 



Comparing the above expression to that of the T' diagram, we have EM, = EM,, 
in the SU(3) limit, leading to the observation that, at 0(a°), the ratio C'/T' is inde- 
pendent of matrix elements [Eq. ()53p]. However, if we combine the 0(q;°) and 0{al) 
results, the ratio is not independent of matrix elements: 



-Nr. 



+ C2) em',, + f (1 - ^) EMa^ 



T' [c, + §-) EM', + f (1 - 1 eLi em 



(64) 
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Note that the 0{al) corrections to C and T' are proportional to ci and C2, respectively. 
Since ci is quite a bit larger than C2, we therefore expect that the first-order correction to 
the ratio C'/T' is also large. However, to obtain the exact value of C'/T' beyond leading 
order, one needs to use a method to estimate the matrix elements. This establishes 
the connection between contractions and QCDfac/pQCD. 



6 Is there a. B ^ ttK puzzle? 



The amplitudes for the four B —>■ ttK decays are given in terms of diagrams in Eq. (j^Uj) . 
Many of these diagrams are expected to be negligible {A', P'j^, P'^w) |2]- Retaining 
only those which are expected to be sizeable, we have 



^+0 


= PL + PLe''- 


-'- p>c 










= -TV^-CV^ 


- P' - 


_ P' p*7 _ 


P' - 

EW 


2 

pic 

tj^^ EW 




= -tV^-p4- 


p' pi'T 


2 

^ pic 








= -cV^ + p;,+ 


PLe'' 


-P' - 

EW 


^ pic 





(65) 



In the above, P4 = P/ — Pc, Pic = PL ~ P'c^ ^'^d we have explicitly written the 
relative weak phase 7 (the phase information in the CKM quark mixing matrix is 
conventionally parametrized in terms of the unitarity triangle, in which the interior 
(CP-violating) angles are known as a, (3 and 7 



The diagrams P^^^ and P'^^ are not independent - as has been shown, they are 



related to T' and C . If we do not make the approximation that Ci/c2 = Cg/cio, 
have 



we 



P' 

EW 



VC 



3 C9 + Clo 

4 Ci + C2 

3 C9 -|- Cio 

4 Ci + C2 



where 



P = 



R{T' + C) 
R{T' + C) 



3 C9 — Cio 

3C9 



- C2 
ClO 



4 Ci - C2 



R{T' - C) 
R{T' - C) 



ub ''us 



1 sin(/3 + 7) 
A2 sin/5 



(66) 



(67) 



(If C1/C2 = C9/C10 is assumed, the ratios of Eq. (jl^ are reproduced.) 

The fact that R can be expressed in terms of /3 and 7 was used in Ref. [H] to 
extract these CP phases from B ttK decays. However, j3 and 7 can also be obtained 
independently: f3 can be taken from the measurement of mixing-induced CP violation 
in P0(t) ^ J/ijKs-. sin 2/3 = 0.687 ± 0.032 while 7 can be found via a fit to 
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Table 2: Branching ratios, direct CP asymmetries Adir, and mixing-induced CP asym- 
metry Amix (if applicable) for the four B ttK decay modes. 



Mode BR{10-^) Adir A 



B+ - 




24.1 ± 1.3 


-0.020 ±0.04 


B+ - 




12.1 ±0.8 


0.04 ±0.04 


B'd- 


^ 7r-K+ 


18.9 ±0.7 


-0.115 ±0.018 


B'd- 




11.5 ± 1.0 


0.02 ±0.13 0.31 ±0.26 



independent measurements: 7 = 58.6l5;9° PH]- In the fits which follow, we include 
these independent determinations of the weak phases. 

Given that P^^^ and P'^y^ can be related to T' and C", the amplitudes in fact 
depend on only seven unknown theoretical parameters: the four magnitudes \Plc\i 
\C'\ and \Puc\i ^^"i their three relative strong phases. However, there are nine B —>■ rrK 
measurements: the branching ratios for the four B — >■ irK decays, the four direct CP 
asymmetries Adir, and the mixing- induced CP asymmetry A^i^ in B^ — > n^K^. The 
latest data is shown in Table El 

It is therefore possible to perform a fit. We find a good fit: Xmin/^-'^-f- = 3.03/2 
[THj . However, the fit also gives \C'/T'\ = 1.73 ± 1.01, whose central value is far larger 
than that given by the result [Eq. This is the B — ^ nK puzzle: present 

B — >■ 7tK data seem to be inconsistent with naive SM predictions. (Above, we used the 
ratio \C'/T'\ to illustrate this point, but other quantities can be used as well.) Many 
analyses have found this result ^7]. Note that, at present, the experimental errors 
are still quite large (and the theoretical errors, such as SU(3) breaking, have not been 
included), so that the effect of the B — >■ ttK puzzle is not yet statistically significant. 
As such, it can be said to offer only a hint of a discrepancy. 

The purpose of this section is to critically re-examine the question of whether there 
is a S — > ttK puzzle. Assuming that the theoretical uncertainties are under control, 
and that the effect is not a statistical fluctuation - and these assumptions may well be 
wrong - the relevant question is: is it possible that \C'/T'\ = 1.7? 

We begin with the expression for C'/T', Eq. (jMj) : 



(t + C2) 


EM',, 


^ 2 






=1 EMai 


+ t 


)EM', 


± ^ 
^ 2 






=1 EMi 



(68) 



The important question is, to what extent does the lowest-order (naive factorization) 
result for C'/T' [Eq. represent the true SM value for this ratio? Note that if 

\C' /T'\ = 0.17 ± 0.1 is included in the flt as a constraint, we obtain a very poor fit: 
Xmin/d-O-f- = 20.6/4, which corresponds to a deviation from expectations of 3.6a (T^ . 
Note also that our fit includes \Puc\y hence a large value of P^ cannot lead to a 
good fit to the data (this is somewhat contrary to Ref. [HI)- 
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In QCDfac the value of \C' /T'\ at 0{a].) may be raised to about 2-3 times the 
lowest-order result jTHj and thus falls far short of the value required to fit the data. 
The fact that the first-order correction is quite large follows from the observation that 
C and T' get corrections proportional to Ci and C2, respectively [Eq. (jMj) ]. but ci is 
about 5-6 times larger than C2. Within pQCD, one also finds that the value of \C' /T'\ 
may be raised to 2-3 times the lowest-order result including NLO corrections PUI. Note 
that although the authors of Ref. ^20j give results that are consistent with the central 
values of the direct CP asymmetries in the Ktc system, they cannot explain the central 
value of the indirect CP asymmetry in K^-k^ . On the other hand, our fit takes into 
account all Ktt data, including the indirect CP asymmetery in K^n^. One therefore 
has to be careful about the claim in Ref [201 that the Kir puzzle is resolved within 
pQCD. 

Although we do not have an all-orders result for \C'/T'\ to compare with the fit, we 
can still make several observations. It is true that higher-order effects, such as those 
from additional gluons, will affect the lowest-order result for \C'/T'\ and change its 
value. However, in order to produce a value that agrees with the fit, this value must 
change by a factor of 10! There are then two possibilities. One is that the true value 
of \C'/T'\ is large, which means that the corrections are very important. In this case, 
there is no 5 — nK puzzle. However, both methods of calculating matrix elements 
(QCDfac and pQCD) have argued for a reasonably small value of the strong couphng 
as at the scale rrib. Thus, the corrections are increasingly small. This in turn 

leads to a perturbative expansion for nonleptonic B decays. Now, if there is really no 
B —>■ ttK puzzle, this means that the 0{a'^) corrections are large, which sheds much 
doubt on the results of both QCDfac and pQCD. The other possibility is that the 
corrections are small and the true value of \C'/T'\ is similar to the 0{a^) result. In 
our opinion, this situation is more likely: we have small corrections, so that the lowest- 
order result for \C'/T'\ [Eq. is approximately correct, the results of QCDfac and 
pQCD are believable, and there is a B ^ nK puzzle. 

7 Conclusions 

We have presented a general approach to hadronic B decays. The starting point is the 
observation that, at the quark level, the decay B —>■ M1M2 (Mi 2 are mesons) involves 
the matrix elements {qiq2qz'ii\'ibb q^qilqsh) , where Mi = qiq2, M2 = ^3^4, and q^hq^qj 
is an operator of the effective Hamiltonian (Dirac and colour structures are omitted 
here). In order to determine all possible ways in which the final-state quarks can be 
produced, it is necessary to sum over all possible Wick contractions of the quarks for 
all operators. The amplitude for B M1M2 can then be expressed as a sum over these 
contractions. In this paper we have examined various properties of these contractions. 

There are a total of 24 possible contractions. However, due to the fact that we can 
write the final state as M1M2 or M2M1, not all contractions are independent. We have 
shown that there are in fact only 14 independent contractions. Buras and Silvestrini 
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(BS) obtained these results some years ago [5]. However, our analysis goes beyond that 
of BS in several ways, described below. 

We have separated the independent contractions into four classes. Using recent 
theoretical and experimental developments, we have shown that certain contractions 
are smaller than others, and can occasionally be neglected. This greatly simplifies the 
expressions for the amplitudes in terms of contractions. 

It is also possible to write the 5-decay amplitudes in terms of diagrams |:_2j. How- 
ever, the relation between diagrams and the effective Hamiltonian was not made clear, 
and there was some question about the rigourousness of the diagrammatic approach. 
We have shown that all diagrams can be simply expressed in terms of contractions, 
thereby demonstrating formally that the diagrammatic method is rigourous. 

In the limit of neglecting the (small) Wilson coefficients C-j and Cg, and assuming 
flavour SU(3) symmetry, Neubert/Rosner and Gronau/Pirjol/Yan have shown that 
there are relations between the electroweak-penguin and tree diagrams. We show that 
these relations are reproduced using the contractions method. 

All of the above results hold to all orders in a^. That is, the contractions method 
includes the exchange of any number of gluons between the quarks. However, we have 
shown that it is also possible to work order-by-order in by including the gluons ex- 
plicitly. This is useful, as it allows one to make a connection between the approach of 
contractions and the matrix-element evaluation methods of QCD factorization (QCD- 
fac) and perturbative QCD (pQCD). If one works to leading order [0(a^)], we have 
shown that one finds that in the SU(3) limit the ratio C jT' is independent of matrix 
elements (T' and C are, respectively, the colour-allowed and colour-suppressed tree 
diagrams in i? — >■ i^K decays). In addition, this ratio is found the be rather small: 
C /T' ^ 0.17. If one adds a single gluon [0(a])], this value changes, but one needs to 
evaluate matrix elements to determine its value. 

Finally, we re-examine the question of whether there is a ttK puzzle." If one 

considers all B — >■ irK data, a good fit is obtained, but \C'/T'\ = 1.73±1.01 is required. 
This central value is far larger than the 0(a°) result. Leaving aside the possibility of 
a statistical fiuctuation, which might indeed be the true explanation, the question of 
whether or not there is a i? ^ ttK puzzle comes down to the question of whether or not 
additional gluons can change \C'/T'\ from 0.17 (the 0{a^) result) to 1.7. The higher- 
order effects have been evaluated in QCDfac and pQCD, and both methods find that 
\C'/T'\ is changed by at most 2-3 times the lowest-order result, instead of the required 
factor of 10. This is not surprising as both methods find that is relatively small at 
the scale rrih. If \C'/T'\ were changed by a large amount, that would shed doubt on the 
smallness of as, as well as all predictions of QCDfac and pQCD. We therefore conclude 
that it is likely that \C'/T'\ is relatively small (i.e. similar to the 0{a^) result), and 
that there really is a i? ^ ttK puzzle. 
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